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Abstract 

The experimental results presented in this paper show the influence of the elastic defor- 
mation of porous solids on the adsorption process. With p"'"-type porous silicon formed 
on highly boron doped (100) Si single crystal, we can make identical porous layers, either 
supported by or detached from the substrate. The pores are perpendicular to the substrate. 
The adsorption isotherms corresponding to these two layers are distinct. In the region pre- 
ceding capillary condensation, the adsorbed amount is lower for the membrane than for the 
supported layer and the hysteresis loop is observed at higher pressure. We attribute this 
phenomenon to different elastic strains undergone by the two layers during the adsorption 
process. For the supported layer, the planes perpendicular to the substrate are constrained 
to have the same interatomic spacing as that of the substrate so that the elastic deformation 
is unilateral, at an atomic scale, and along the pore axis. When the substrate is removed, 
tridimensional deformations occur and the porous system can find a new configuration for 
the solid atoms which decreases the free energy of the system adsorbate-solid. This results 
in a decrease of the adsorbed amount and in an increase of the condensation pressure. The 
isotherms for the supported porous layers shift toward that of the membrane when the layer 
thickness is increased from 30 to 100 ^m. This is due to the relaxation of the stress exerted 
by the substrate as a result of the breaking of Si-Si bonds at the interface between the 
substrate and the porous layer. The membrane is the relaxed state of the supported layer. 

Author E-mail address : annie.grosman@insp.jussieu.fr 

1 Introduction 

The elastic deformation of porous materials during the gas adsorption-desorption process has 
been studied for a long time. 1 The reversible formation of an adsorbed film preceding capillary 
condensation leads to an extension of the porous material, while, when condensation occurs, 
a contraction is observed, attributed to the negative fluid pressure under the concave menisci. 
As the vapor pressure is further increased up to the saturation vapor pressure, an extension 
occurs due to the vanishing of the negative liquid pressure consequent upon the flattening of the 
concave menisci. During the evaporation process, the hysteresis phenomenon is observed in the 
same pressure region as the adsorption phenomenon. [51 13 HI [3 [S] Similar results have been 
obtained in SBA-15 ^ and in MCM-41 ^ by means of in situ small-angle X-ray diffraction 
experiments. 

Curiously, these deformations were generally ignored in the thermodynamics of adsorption. 
How does the elastic deformation of a porous solid affect the adsorption and capillary condensa- 
tion? In the recent paper quoted above, [S] the authors address the question by means of Monte 
Carlo simulations. They found that the elastic deformation causes a shift of the order of a few 
percent of the phase transition toward that in the bulk. 

Simultaneously with this work, [8] using a classical thermodynamic approach, we have shown 
[2 that the assumption of an inert adsorbent in adsorption is totally unjustified. In the physics 
of adsorption, the equilibrium condition for mass transfer is not given by the equality of the 
chemical potentials of the vapor and adsorbate but includes also the deformation of the solid. 
This allows us to introduce in the thermodynamic relationships, in addition to the classical surface 
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free energy, the elastic energy stored in the solid, which is, for small deformations, quadratic in 
the deviation. The mechanical equilibrium obtained by the minimization of the free energy of 
the adsorbate-solid system with respect to the deformation results in the linear variation of the 
deformation with the surface free energy observed by some investigators, [3l I10|, lllj with the 
proportionality factor depending on the elastic constants of the solid. 

The linear variation of the deformation with the surface free energy is hence the result of two 
equilibrium conditions, with one being thermodynamic and the other mechanic. This important 
result means that the adsorbed amount in a given pore can be changed by a mechanical stress 
external to the pore. 

This is the basis of the interaction mechanism we have proposed. The filling or the draining 
of a pore causes elastic deformation not only on its inner pore walls but also on the inner walls of 
its neighbors and thus causes a change of their surface free energy, i.e., a change of the adsorbed 
amount in a neighbor pore not yet filled and a change of the negative liquid pressure through 
the shape of the meniscus in a neighbor pore already filled. This interaction mechanism acts 
as an alternative to the pore blocking/percolation mechanism, [12l [13] which fails to explain 
the interpendence between pores observed in systems where there is no pore-pore intersection. 

[Hill] 

In the present paper, we find the same idea as that exposed in our interaction mechanism 
except that here the external mechanical stress is not imposed on a given pore by the neighbors 
but imposed by the substrate to which the porous layer is attached, which is the case for porous 
silicon. We thus come back to a problem which we have already addressed in a letter: 
the comparison of the adsorption isotherms for p"*"-type porous silicon layers supported by the 
substrate with pores closed at one end and detached from the substrate (a membrane) with 
pores open at both ends. Capillary condensation systematically occurs at higher pressure in the 
membrane than in the supported layer. Following the idea of Cohan, [l?] we attributed this 
phenomenon to two different scenarios for the filling of these two porous systems. In a pore 
closed at one end, a preferential condensation occurs at the bottom of the pore, while in a pore 
open at both ends the adsorbed film becomes unstable and the filling occurs at higher pressure 
when a liquid bridge is formed somewhere in the pore. 

In the light of our thermodynamic approach [3] and of new experimental results presented 
in this paper, we reconsider this explanation. After discussing the similarity of the supported 
layer and the corresponding membrane in Section [3.H we show in Section [3. 21 that the difference 
between the two isotherms is due to the different elastic deformations undergone by the two 
layers on account of the stress exerted by the substrate: for the supported layer, the deformation 
is unilateral (perpendicular to the substrate) while the membrane is submitted to tridimensional 
deformations not necessarily isotropic. To these different configurations for the solid atoms 
correspond different free energies. In Section [3. 3 [ we show that this effect disappears gradually 
when increasing the layer thickness from 30 to 100 /im: the isotherm corresponding to the 
supported layer shifts toward that of the membrane. This is due to the relaxation of the stress 
exerted by the substrate as some Si-Si bonds are broken at the porous layer-substrate interface 
by the stresses generated at this interface. In Section [4] we will discuss these whole results. 

2 Experimental Section 

2.1 Preparation and Properties of p+-type Porous Silicon 

The preparation process, the morphology, and the physical properties of p+-type porous silicon 
have been reported in detail in recent papers. [lH [15] We give here a brief outline of these 
properties. p"'"-Type porous Si is obtained by anodic dissolution, at constant current, of highly 



2 



boron doped (100) Si single crystal in hydrofluoric acid (HF) solution. It is a porous single 
crystal material. The pores are perpendicular to the Si substrate, and all have the same length. 
Transmission electron microscopy (TEM) measurements show that (i)- at the surface of the 
porous layer, the pore cross sections are of polygonal shape, (ii)- the pore size distribution 
(PSD), which depends on the formation conditions (current density and HF concentration), is 
large 13 ± 6 nm and ^ 26 ± 14 nm for a 50% and 70% porosity layer, respectively), and 
(iii)- the thickness of the Si walls separating the pores is almost constant. The two-dimensional 
(2D) TEM images yield the same porosity as that obtained by weighing. Moreover, the pore 
volume measured by weighing or by adsorption measurement once the pores are full of liquid is 
proportional to the etching time. This indicates that during the formation process the Si walls 
become insulating once a certain thickness is reached, 5 — 6 nm in our case, even though the 
dopant atoms are not removed from the porous layer during the formation process. [TH] Thus, 
except for the presence of numerous facets on the Si walls, [H] i.e., a roughness at the scale of 
few atoms, the pore cross-sectional areas are constant in depth. The cause of this high resistivity 
is not yet well understood. What we know is that the Si-HF interface plays a fundamental role 
in this property, since the Si walls remain electrically resistive as long as the HF environment is 
present. To render the Si walls electrically conducting, we must replace the acid HF solution by 
a neutral solution; thus, anodic oxidation occurs in place of anodic dissolution. This electrical 
property can be used to make membranes. At the end of the preparation process, once the desired 
pore length is obtained, the anodic current can be either cut off to form a porous layer supported 
by the Si substrate with pores closed at one end or increased to go into an electropolishing regime 
during which the Si walls are dissolved at the bottom of the pores to form a membrane with pores 
open at both ends. For this paper, we have prepared porous layers with a porosity of 50% and of 
different thicknesses varying between 10 and 100 /xm. We have also prepared the corresponding 
membranes. The electrochemical conditions were HF(40%)/EtOH=3:l, J = 20 mA/cm^, and 
etching rate ~ 1 /Ltm/min. The main structural parameters of this kind of layer can be found in 
ref 15. 

2.2 Measurements 

The nitrogen adsorption isotherm measurements were made at 77.4 K using a Micromeritics 
ASAP2010 instrument equipped with pressure transducers with full scales of 1 /imHg, 10 mniHg, 
and 1000 mmHg. 

The radius of curvature of the supported porous layers was determined optically by measuring 
the diameters of Newton's rings. For these measurements, we have prepared porous layers on 
Si wafers polished on both sides. The interference patterns were caused by the difference in 
path lengths of the reflecting light from the concave back side of the Si substrate and from a 
flat mirror, with the whole device being in the configuration of a Michelson interferometer with 
equal arms. We used a red laser light (wavelength of 0.6328 //m). 

The cross section view of the porous layer-silicon substrate interface was obtained using the 
high spatial resolution of a Field Emission Gun Scanning Electron Microscope (FEG-SEM). The 
beam energy was 5 kV and the magnification was 200000 x . 
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Figure 1: Nitrogen adsorption isotherms at 77.4 K for several p'''-type porous silicon layers 
supported by the Si substrate (symbols: •, ♦) and the corresponding membranes (symbols: 
□ , o, A, 0)- The solid lines are interpolations of the mean values. The thickness equals 20 /im, 
and the porosity is 50%. 

3 Results and Data Analysis 

3.1 The porous layers supported by the substrate and the correspond- 
ing membranes 

As noted in the Introduction, we come back, in this paper, to a problem which we have 
already addressed: the comparison of the adsorption isotherms for porous layers supported by 
the substrate with pores closed at one end and for self-supported layers with pores open at 
both ends. Figure [T] shows nitrogen adsorption isotherms at 77.4 K, normalized to the same 
pore volume, for several 20 /im thick supported layers and the corresponding membranes, with 
porosity of 50%. The solid lines are interpolations of the mean values of the two types of layers. 
Note that when making membranes, it is difficult to collect all the sample pieces, with some of 
them being too small and others, located at the periphery of the layer, remaining attached to 
the substrate. It is the reason why the isotherms for the supported layer and the corresponding 
membrane were normalized to the same pore volume. 

The condensation branch was systematically observed at higher pressure for the membrane 
than for the supported porous layer. Following the idea of Cohan, [T7] we attributed \l6f this 
phenomenon to the effect of pore ends on the condensation process. 

However, two results shown in Figure [T] remained unexplained. (1) In the region of reversible 
adsorption preceding the capillary condensation (p ^ 0.5), the adsorbed amount is systematically 
lower for the membrane than for the supported layer. (2) The evaporation branches of the layer 
and the corresponding membrane are not superimposed. 

Under the usual assumption of an inert adsorbent, the presence of a second opening in 
otherwise identical porous layers, where the pore cross section is invariant with depth, should 
change neither the adsorbed amount before condensation nor the position of the evaporation 
branch on the pressure axis. 
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Figure 2: Nitrogen adsorption isotherms at 77.4 K for 100 /im thick supported porous sihcon 
layers (symbols:*, ▲) and the corresponding membranes (symbols:o, □). The dashed and solid 
lines are, respectively, the interpolations of the mean values for the 20 /im thick supported porous 
layers and membranes represented in Figure [T] The porosity is 50%. 



As the results of Figure [T] are the central point of the paper, we find it important to come 
back to the problem of the strict similarity of the supported and detached layers. 

The only difference between these two layers, if any, could be a shift of the whole PSD 
equal to a fraction of the pore wall thickness, for example, ~ 1 nm, under the effect of the 
electropolishing process. This process, which dissolves the Si walls at the bottom of the pores, 
could slightly dissolve the Si walls in the whole porous system, increasing the pore size. The 
determination of the PSD using TEM, described in detail in refs 15 and 16, is not sufficiently 
precise to measure such a shift. Our PSD measurements concern the supported layers. 

We have shown [T5j that the mean pore diameters defined by < Dp >=< p > /it, where 
< p > is the mean pore perimeter, and by < Dg >= (4s/7r)^/^, where < s > is the mean surface 
area of the pore section, are quite similar (< Dp >= 13.2 nm and < Ds >— 13.3 nm). Thus, 
the S/V ratio of the pore surface area S — NL < p > over the pore volume V = NL < s > is 
inversely proportional to the common value < D > (N is the number of pores, and L is the pore 
length) . The adsorbed amount before capillary condensation normalized to the same pore volume 
is hence inversely proportional to < Z? >, so that, if we suppose that the deformation of the 
adsorbent does not play any role on the adsorption process, the comparison of the normalized 
adsorbed amounts before capillary condensation for the supported layer and membrane is an 
excellent criterion to decide whether these two layers are identical. The lower is the normalized 
adsorbed amount before capillary condensation, the higher is the mean pore size. 

Thus, under the assumption of an inert adsorbent, the conclusion of the results of Figure [T] is 
that the mean pore size of the 20 fun thick membrane is higher than that of the corresponding 
supported layer. Concerning the condensation process, the shift of the condensation branch 
toward higher pressures observed for the membrane could be explained either by the sole increase 
of the pore size without invoking any influence of pore end on the condensation process or by 
the addition of the two effects. 

Nevertheless, the above explanation seemed to us very unlikely. The electrochemical etching 
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rate is typically a few /im/min and the Si walls are 5 — 6 nm thick. The simple fact that we 
can obtain membranes of different thicknesses by this process, with reproducible normalized 
adsorption-desorption isotherms (see below), without paying special attention to the value of the 
current intensity chosen in the large electropolishing regime, ^ 100 mA/cm^ in our experimental 
conditions, strongly suggests that the porous layer remains unchanged during the electropolishing 
process. 

Fortunately new results made it possible to solve this problem. Indeed, when making thicker 
supported layers (100 /im) and the corresponding membranes, we fortuitously fell upon a surpris- 
ing result. As shown in Figure [21 the adsorption-desorption isotherm corresponding to a 100 /xm 
thick supported layer is fully identical to that of the corresponding membrane, and in addition 
to that of the 20 fim thick membrane. The experiment was repeated several times with similar 
results. 

Two conclusions can be drawn. 

(i) The isotherms corresponding to the thin and thick membranes are superimposed once 
normalized to the same pore volume, which shows that they differ only by their thickness. If the 
electropolishing process used to obtain the two membranes modifies the PSDs, the modification 
is independent of the porous layer thickness. In the region preceding capillary condensation, the 
adsorption isotherms corresponding to the 100 fim thick supported layers are superimposed on 
those corresponding to the two membranes. Making the same assumption of an inert adsorbent 
as for Figure [1] we conclude that the electropolishing process lets the 100 /im thick supported 
layers unchanged. This is hence also the case for the 20 /im thick supported layer so that it is 
similar to its corresponding membrane. This is in contradiction with the above interpretation of 
Figure [TJ 

(ii) The isotherms are also superimposed in the region of capillary condensation, which means 
that the thick supported layers behave like the membranes, which is not the case for the thin 
supported layers. Concerning this point, the question is whether this could be simply explained by 
the disappearance of pore end effect on the condensation process when making thicker supported 
layers. If we look at the change of the adsorption isotherm when increasing the thickness of the 
supported layer from 20 to 100 /im we note that the shift of the condensation branch toward 
higher pressures is accompanied by a decrease of the normalized adsorbed amount before capillary 
condensation. This indicates that this shift is not only due to pore end effect, if any. Indeed, 
still under the assumption of an inert adsorbent, pore end effect can influence the condensation 
process but should not affect the adsorbed amount before condensation. 

We conclude that the assumption of an inert adsorbent during the adsorption process must 
be abandoned if we want to explain the whole results shown in Figures [T] and [51 So, in the 
following section, we will analyze the above striking results by taking into account the elastic 
deformation of the porous matrix. At this stage, since the adsorbent cannot be regarded as 
inert anymore, the comparison of the different isotherms does not allow us to draw a conclusion 
concerning the similarity of the supported layers and the corresponding membranes. In what 
follows, we have chosen to assume their strict similarity, a hypothesis which seemed to us by far 
the most plausible and which will be validated a posteriori. 

In Section 13.21 we first analyse the results shown in Figure [H and, in Section 13.31 those of 
Figured The possible effect of pore ends on the condensation process will be discussed. 

3.2 Elastic deformation of the porous layers 

The lattice parameter of porous silicon layers is different from that of silicon: [50] porosity 
affects the crystal lattice by producing a slight expansion. For p+-type porous layers similar to 
those we study (highly boron doped (100) Si wafer, porosity ~ 50%), the lattice parameter of 
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porous layers supported by the substrate (subscript si), measured in the [100] direction, is bigger 
than that of the substrate: {Aa/a)si = 5.15 x 10"'*. For the corresponding membrane (subscript 
to) the lattice parameter (still in the [100] direction) is a little less large: (Aa/a),„ ~ 4.2 x 10"'*. 
This is due to the fact that, for the supported porous layers, the planes perpendicular to the 
substrate are constrained to have the same interatomic spacing as that of the substrate, leading 
to an additional expansion: for the layers studied by Barla ct al., [50] {Aa/a)si — 1.2 {Aa/a)m- 

The expansion of porous silicon is attributed to the presence of SiH2:=i^2,3 groups chemically 
adsorbed during the formation process, [211 [12] but, to our knowledge, nobody has measured the 
lattice parameter of clean porous layers. Once the layers are exposed to air, the hydrogen atoms 
are slowly replaced by oxygen over a timescale of a few months. [23] The two (Aa/a) values 
given above probably correspond to freshly prepared porous layers covered by the same amount 
of SiHa; groups, i.e., to the same surface stress. 

Thus, during adsorption experiments, for a given adsorbed amount, we should have 

-pC-^) , (1) 



where Sa is the variation of the lattice parameter (in the [100] direction) caused by adsorption 
and /3 > 1 is a constant which depends on the morphology of the porous layer. /3 = 1.2 for the 
layers of ref 20. 

We have no information on the elastic deformation {5a/ a) of the supported porous layer 
and the corresponding membrane during N2 adsorption experiment, but several investigators 
[Sidnilll] have shown that, in the region preceding capillary condensation, the linear deformation 
of porous materials is related to the adsorbed amount by the equation 

N^kBT'^^k-, (2) 
P a 

where N„ is the number of the adsorbed molecules, P is the vapor pressure, and fc is a propor- 
tionality factor which depends on the elastic constants of the porous solid. 
We have shown ^ that 
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, p dP 
lA\l = -j^ N^ksT—, (3) 



where is the elastic energy of the solid, 7 is the surface free energy per unit area, and A is 
the surface area of the porous solid. The elastic energy stored in the solid during adsorption can 
be neglected compared to the variation of the surface free energy so that eq 3 becomes 

,p dP 
A{jA) = \jA\^ = - N^ksT— (4) 

which is the classical Gibbs adsorption equation established under the assumption of an inert 
adsorbent. From eqs 3 and 4, it follows 

\jA\^ =[j{a)-j{a,)]A = -k-, (5) 

I lU L 

where ag and a are the lattice parameters before and after adsorption, respectively. 

The variations of the surface free energy as a function of the vapor pressure, corresponding 
to the two isotherms represented in Figure (TJ in the region preceding capillary condensation, are 
shown in Figure[3| For a given adsorbed amount, i.e., a given variation of the surface free energy. 
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Figure 3: Variation of the surface free energy as a function of the gas pressure, deduced from ni- 
trogen adsorption isotherms at 77.4 K in the region preceding the capillary condensation process, 
for a 20 /xm thick supported porous silicon layer (thick line) and its corresponding membrane 
(thin line). The mass of the porous layer equals 36.36 x 10"'^ g. 



the elastic deformations of the supported layer and membrane are related by eq 1. It follows 
from eq 5 that 

km = l3 ksu (6) 

where km and ksi are the proportionality factors, defined by eq 5, corresponding to the membrane 
and the supported layer, respectively. 

The ratio A (73; A) /A (7mA) calculated from the data represented in Figure [3] is independent 
of the vapor pressure and equals 1.111 ± 0.001. 

Thus, at each vapor pressure during the adsorption process, the extension of the supported 
layer and the membrane should obey the relation 

Taking /3 = 1.2, for example, we obtain {5a/a)si = 1.33 {5a/a)m- A decrease of ~ 30% in {Sa/a) 
along the pore axis accompanied by transversal deformations then leads to a significant difference 
(10%) in the adsorbed amounts (see Figure[T]). 

It remains now to verify whether the variation of the surface free energy for either porous layer 
corresponds to acceptable values of the relative deformation. It is easier to make this estimation 
for the supported layer, since the deformation is unilateral. We have shown [5] that k — CVso/a, 
where C is defined by 

C is a constant which depends on Young's modulus Ep, on Poisson's coefficient Vp, and on the 
geometry of the porous solid. Vso is the porous solid volume, and a is a constant defined by 
dA/A = a{da/a). 



8 



Equation (5) becomes 

[7(a)-7M]^ = -^ (-)■ (9) 
a \ a J 

For the supported porous layer, perpendicular to the pore axis {Sa/a)± = 0. In this case, 
dA/A = da /a and a = 1. According to the classical theory of elasticity, the constant C is 
related to Young's modulus by the relation 

As ~ 0.1, [30] C ~ £'p and hi ~ Ep V^o. 

Young's modulus values Ep of p"'"-type porous silicon measured using different techniques |25j 
are close to E' (1 — PorY , a relation found by Gibson and Ashby, [26] where E = 166 GPa is 
Young's modulus of bulk silicon [57] and Por = 0.5 is the porosity of our porous layers. Taking 
a density of 2.32 g/cm^ for Si, we find ksi = 1.8 x 10"* J/g. Considering the data represented in 
Figured we have, at P = 5x10* Pa for instance, A(7^;A) = -0.21 J/36.36 x 10"^ 5 = -5.78 J/g, 
which corresponds to {6a/ a) si — 3.2 x lO"'^. 

It must be noted that the elastic deformations measured before or during adsorption are 
caused by the interaction adsorbate-porous solid and must be considered at an atomic scale. 
Hence, in adsorption experiments. Young's modulus may be higher than that given by the Gibson- 
Ashby relation which corresponds to the elastic response of a porous solid to an external loading, 
and the actual deformations could be lower than those estimated in this paper. 

We now analyze the above results. 

(1) First of all, let us compare the surface free energy to the other terms of the total free 
energy of the system adsorbate-porous solid which is given by [5] : 

Fa^so = -PaVa + H^N, + Ff^ + 7A. (11) 

where ^A ^ F^^ depends on the elastic deformation through the mechanical equilibrium eq 5. As 
long as we are concerned by the reversible adsorption of a few monolayers the term Po-Vo-, where 
Po- and Vc are the pressure and volume of the adsorbed phase, is much smaller that the two other 

terms and can be neglected. For the supported layer, \Na /Xv|p^^° = \NrjkBT /n(P)|g^^" = 
1.21 J which shows that the corresponding variation of the surface free energy, —0.21 J, cannot be 

I 1 5 X 10^ I 1 5 X 10^ 

neglected. Concerning the membrane, \NakBTln{P)\^ = 1.08 J and |7^|g = —0.19 J. 
We see that the variation of the total free energy, Po-.so, of the system adsorbate-porous solid in 
the pressure range — 5 x lO'' Pa is lower for the membrane, 0.89 J, than for the supported layer, 
1 J. This can be interpreted as follows: in its relaxed state (the membrane) the porous system 
has three degrees of freedom and can find a new configuration for the solid atoms which decreases 
the total free energy (~ 10%) and the adsorbed amount (10%) and increases the condensation 
pressure with regard to the supported layer. 

(2) We now compare the estimated value of the coefficient ksi to other experimental results 
found in the literature. For adsorption of water in porous glass, [31 [3] the variation of the 
surface free energy in the region preceding capillary condensation equals ~ —17 J/g and the 
corresponding relative extension, measured by Amberg et ai, [3] is Al/l ~ 1.2 x 10~^ which 
leads to fc = 1.4 X 10"* J/g. The k values estimated in a porous silicon layer from theoretical 
considerations independently of adsorption data, kgi = 1.8 x 10^ J/g, and in porous glass from 
the adsorption data and relative extension measurements, k — 1.4 x 10^ J/g, are hence similar. 
This indicates that for both materials the relative deformations during adsorption should be 
approximatively in the same ratio as the surface free energy variations. Note that if porous glass 
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is assummed to be submitted to isotropic strain the coefficient k is related to Young's modulus 
of porous glass Epgi and to Poisson coefficient Vpgi by the classical relation k = (3/2) [Epgi/{1 — 
'^h'pgi)] Vso- Inserting Vpgi ~ 0.2 [21] we find Epgi = 14 GPa while for porous silicon Ashby 
relation gives Ep ~ 41 GPa which is coherent if we consider that bulk silicon is stiffer than bulk 
glass. [251 

For pentanc in porous silicon, [3] the relative extension, 5a/ a, is small, a few 10^'^, before 
condensation occurs, while in the hysteresis region it equals a few 10^**. This experimental 
result is qualitatively similar to that obtained for water in Charcoal [2] : in both cases the elastic 
deformations arc much more important in the hysteresis region. For water in Charcoal this 
behavior is clearly related to the shape of the hysteresis: the adsorbed amount before capillary 
condensation is very small. This should also be the case for pentane in porous silicon. In our 
experiments, the amount of nitrogen molecules adsorbed before condensation nearly corresponds 
to half of the total pore volume (see Figured]). This can explain why the value calculated above, 
{5a/a)si = 3.2 x 10"'*, is higher than that found for pentane before condensation. 

(3) Let us comment on the relation km = P ksi and the value (3 — \.2 found by Barla et al. 
|20j . For supported layers we have seen that kgi ~ Ep Vso- For membranes, under the assumption 
of isotropic deformation, we should have fc„j = (3/2) [Ep/{1 — 2i>p)] Vso- By inserting Vp ~ 0.1 
we find k,n ~ 1.9 kgi- Note that the low value i/p = 0.1 corresponds to a solid submitted to 
high anisotropic strains. For isotropic material Vp is certainly higher so that the actual value 
of P may be even larger than 1.9. Similar f3 values have been found for cpitaxic layer on (100) 
substrate for which the lattice mismatch measured in the [100] direction has the same value in 
the two other directions. For example, for epitaxic layer Al^Gai_2;As on (100) GaAs substrate, 
for a lattice mismatch of 4.86 x lO"'', /3 = 2.2. [53] 

The fact that, for porous silicon, the coefficient (3 is close to 1 indicates that this porous 
material does not behave elastically neither as a cubic epitaxial layer on silicon nor as an isotropic 
porous medium. It is rather an orthotropic system in which the properties are the same in the 
(Ox) and (Oy) directions parallel to the interface but different in the (Oz) direction perpendicular 
to the interface. The deformations perpendicular to the pore axis (transversal deformations) are 
much smaller than the deformations parallel to the pore axis (longitudinal deformations). 

Let us note that, in ref 20, the Poisson's coefficient value i^p ~ 0.1 has been deduced from the 
relation [30] 

Aa\ ^ /AaX ^1+^/AaX ^^^^ 

This relation, based on the assumption of a homogeneous material, cannot be applied to p"'"-type 
porous silicon. Nevertheless, if we consider the high anisotropy of the material, such a small 
value around 0.1 seems quite plausible. 

We can give a qualitative description of the deformation of the membrane compared to 
that of the supported layer. Parallel to the pore axis (Oz) the extension is probably uniform 
for both layers. In the plane perpendicular to the pore axis, the interatomic spacing in the 
supported porous layer is the same as that of the substrate while for the membrane the situation 
is more complicated. In this plane, for a given pore wall, two directions must be considered 
(see Figure H^). Parallel to the pore walls (Ox), the surface stress exerted by the chemisorbed 
molecules (e.g., SiH^; groups) increases the interatomic spacing and thus the pore perimeter with 
regard to those of the supported porous layer. Perpendicular to the pore walls (Oy), the lattice 
strains are different. Indeed, the contraction of the pore walls along (Oz), compared to the 
supported layer, involves an extension of the pore wall along {Oy), proportional to Poisson's 
coefficient. On the other hand, the extension of the pore wall along (Ox) involves a contraction 
along (Oy). Considering that the contraction of the pore length and the increase of the pore 
perimeter have opposite effects on the pore wall thickness and that Poisson's coefficient is small. 
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(a) (b) 




Si 



Figure 4: Schematic representation of the effects of the relaxation of the stress exerted by the 
substrate on a pore in a plane perpendicular to the pore axis when detaching a supported porous 
layer (a) to obtain a membrane (b): the pore perimeter increases, and the wall thickness is 
unchanged. The lattice mismatch between the membrane and silicon produces an elastic bending 
of the substrate (c), which becomes convex with a radius of curvature R. 
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we can conclude that the wall thickness should be similar for both layers. The effects of the 
relaxation of the stress exerted by the substrate on a pore are schematically represented in 
Figures UK and b. In summary, perpendicular to the pore axis, the elastic deformation in the 
relaxed state of porous silicon (the membrane) is not uniform. Removing the substrate increases 
the pore size but lets the pore wall thickness remain unchanged. It is perhaps why the (3 value 
is close to 1. 

As seen below, the stresses at the interface porous layer-substrate produce an elastic bending 
of the substrate which becomes convex (see Figure Ht) . The supported porous layer behaves like 
an epitaxial layer over silicon which would have a positive lattice mismatch. We shall come back 
to this problem in the next section. 

Let us conclude this section now. 

(1) The supported porous layer and the membrane react differently to the spreading pressure 
exerted by the adsorbed molecules. The deformation is unilateral (along the pore axis) for the 
supported layer, while the membrane is also submitted to transversal deformations. To this new 
configuration for the solid atoms corresponds a lower total free energy which results in a decrease 
of the adsorbed amount for the membrane compared to the supported layer. 

(2) One can explain why the condensation pressure is higher for the membrane than for the 
supported layer without assuming any effect of a pore end on capillary condensation. Let us 
assume that the pores of the supported layer fill in the same way as the membrane. As the 
vapor pressure is increased, the adsorbed film reaches a critical thickness above which it becomes 
unstable. A liquid bridge is then formed somewhere in the pores, leading to the filling of the 
pores; however, since, for a given vapor pressure, the adsorbed film thickness is lower for the 
membrane than for the supported layer, the filling of the pores will occur at higher pressure, as it 
is observed. The experimental results presented in the next section support this idea according 
to which the gap between the two condensation branches is only due to the different elastic 
deformations to which the supported layer and membrane are submitted. 

(3) Consider the two evaporation branches in Figure [1] they are not superimposed, whereas 
the two layers are identical and the pore cross section does not vary with depth. As already noted, 
under the assumption of an inert solid and whatever the evaporation mechanism proposed, the 
presence of a second end should not change the position of the evaporation branch. The fact 
that the supported layer and the membrane empty at different pressures is hence due to the 
different elastic deformations undergone by the two layers. This supports our thermodynamic 
approach [3] according to which the energy barrier to evaporation is essentially due to the elastic 
deformation of the porous matrix. 

Recently, Bruschi et al. |31] studied gas adsorption in nanoporous alumina which is charater- 
ized by a regular arrangement of noninterconnected cylindrical pores. For 25 nm pore diameter, 
the authors report on a similar result as that shown in Figure [U i.e., the hysteresis loop ob- 
tained in pores open at both ends is shifted to a higher relative pressure than that observed for 
closed-bottom pores. They attribute this shift to a slight widening of the pores when preparing 
the porous membrane. By considering the present study, their results could be explained in the 
same way as ours. 

3.3 Stress at the porous layer-compact silicon interface 

The fact that the 100 /im thick supported porous layer behaves like a membrane during 
the adsorption-desorption process (see Figure ^ suggests that the thick layers are imperfectly 
supported by the substrate, that is, Si-Si bonds break at the bottom of the pores under the 
effect of stresses generated at the interface between the porous layer and the substrate, leading 
to the formation of dislocations. However, as shown below, the layer is not detached from the 
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Figure 5: Nitrogen adsorption isotherms at 77.4 K for 100 /ini thick porous sihcon layer 
(symbol:*) and a 100 /xm thick porous silicon layer chemically etched down to ^ 20 /im 
(symbol: a), both supported by a silicon substrate. The lines are guides for the eye. 

substrate. 

First of all, we have verified that the breaking of Si-Si bonds occurs during the formation 
process and not during the adsorption experiments under the effect of additional stresses or of 
low temperature (77 K), and so on. We have made a 100 ^m porous layer and thinned it in a 
NaOH/H20 solution to obtain a 20 /im thick layer. Let us recall that dried porous silicon is 
hydrophobic so that the solution does not enter the pores, which allows us to thin the porous layer 
without changing the pore size. Figure [5] shows the corresponding isotherm together with that of 
a 100 /im supported porous layer. Once normalized to the same pore volume, the two isotherms 
are superimposed, which shows that the Si-Si bonds break during the formation process. The 
formation of thick porous layers (e.g., 100 /im) changes irreversibly the porous layer-substrate 
interface. 

For what physical reasons do Si-Si bonds break at the porous silicon-substrate interface when 
making thick porous layers? The formation of a porous silicon layer in a silicon wafer produces an 
elastic bending of the sample, which becomes convex. [20] As already noted, the system behaves 
like an epitaxial layer over silicon which would have a positive lattice mismatch, (Aa/a)m±, 
perpendicular to the pore axis. The substrate radius of curvature is given by ^32j 

R= {d + dp) 

where 2d = t and 2dp — tp are the thicknesses of the substrate and of the supported porous 
layer, s — [1 ~ v)/E and Sp = (1 — Vp)/Ep where E,Ep and Vji^p are Young's moduli and 
Poisson's coefficients of silicon and of the porous layer, respectively. The porous layer thickness 
is measured by weighing as described in ref 15. Remember that the elastic deformation in the 
relaxed state of porous silicon (the membrane) is not uniform so that (Aa/a)mi_ is not really 
a lattice mismatch as is the case for epitaxial layers, but only a phenomcnological parameter: 
{Aa/a)m± can be considered as an "effective" lattice mismatch. 

In our case, we have (t + tp)^ 300 /xm and tp < 100 /im. In these conditions, as verified 
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below, eq 13 can be replaced with a good approximation by 
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The difference ((72 — cri ) between the surface stresses exerted on the front side of the wafer where 
the porous layer is formed and on the back side is given from classical elasticity theory by 

6 [I-V) ((72 - (7l) 



From eqs 14 and 15, we get 
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((72 — (7i) is hence proportional to the thickness of the porous layer. Increasing the porous 
layer thickness increases the stresses at the porous layer-substrate interface and can cause Si-Si 
bonds to break. This effect is illustrated in Figure El which shows N2 adsorption isotherms for 
porous layers of increasing thickness from 10 to 100 /im in the hysteresis region (Figure [6^) 
and in the low pressure range (Figure [Hb). Up to a thickness of ~ 30 /im, the whole isotherms 
are superimposed, then they shift toward that corresponding to the membranes, the adsorbed 
amount in the region of reversible adsorption preceding the capillary condensation decreasing 
continuously (see Figure Eb). These results indicate that Si-Si bonds begin to break only when 
the thickness exceeds a limit between ~ 30 /im and ~ 40 /im. From this limit, as the thickness 
is increased, more and more Si-Si bonds break so that the stress exerted by the substrate on 
the porous layer is less and less observable in adsorption experiments. Consequently, the broken 
Si-Si bonds are not localized only at the porous layer-substrate interface but also in the whole 
porous layer from a depth > 30 /im. 

To confirm this explanation, we have determined optically the radius of curvature of the 
different samples by measuring the diameters of Newton's rings. The results are shown in Figure[7] 
as a function of the porous layer thickness. We have found that the Newton's rings are circular 
for thick porous layers (50 and 100 /im) but elliptical with two orthogonal principal axes for 
thinner layers. For example, for the 20 /im thick layer, we found Ri 3.5 m and i?2 ^ 4 m. The 
R values in Figure [7] correspond to several measurements with a deviation of ~ 2%. As expected, 
the radius of curvature of porous layers decreases with increasing thickness. 

We have checked whether these measurements can be fitted by the theoretical predictions. 
Since Si-Si bonds do not break in the 20 and 30 fim thick porous layers, the curves represented 
by eqs 13 and 14 must pass by the points corresponding to these two layers. In the case of 
eq 14, the only parameter which we need, {sp/s){Aa/a)^^, has been deduced from the radius 
of curvature of one or the other of the two layers. In the case of eq 13, we need to know 
separately the values of the two parameters, (sp/s) and (Aa/a)~^. As we have seen above, 
porous silicon is an anisotropic material. In eqs 13 and 14, the elastic compliance parameter (sp) 
is related to deformations parallel to the interface. This parameter which can be different from 
that corresponding to deformations perpendicular to the interface is unknown. However, we have 
estimated (sp) by assuming that Ep = E (1 — Por)^, Vp ~ 0.1 which corresponds to deformations 
perpendicular to the interface. For bulk silicon, v — 0.26. [20] The second parameter, (Aa/a)„ix, 
has been deduced as above from the radius of curvature of the 20 /xm thick layer. We have found 
[Aa/ a)mi_ — 6.5 x 10^^. If our qualitative description represented schematically in Figure [His 
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Figure 6: Nitrogen adsorption isotherms at 77.4 K for 10 //m, 20 /im, 30 /im, 40 /im, 50 /im 
and 100 thick supported porous layers in (a) the hysteresis region and (b) the pressure range 
[0 — 0.5]. The isotherms for layers 10, 20, and 30 /^m (symbols:*) thick are superimposed in the 
whole pressure range. The lines are guides for the eye. 



15 




Figure 7: Radius of curvature R measured from Newton's rings (symbol: A) and calculated from 
eqs 13 (symbol; o) and 14 (symbol:*) as a function of the porous layer thickness, tp. The R values 
correspond to several measurements with a deviation of ^ 2%. The inset shows the experimental 
values of t^/R as a function of tp, where t is the substrate thickness (symbol: A). The straight 
line which passes by the origin and the first two experimental points is calculated from eq 14. 
The porosity is 50%. 



good, this value should represent an estimate of the extension mean value of the pore diameter 
when the stress exerted by the substrate is relaxed. 

The two theoretical curves corresponding to eqs 13 and 14 are represented in Figure [7] As 
noted above, we see that eq 14 is a good approximation of eq 13. The radius of curvature of layers 
thicker than 30 fiui deviates more and more, in relative terms, from the theoretical values, which 
shows that an increasingly large number of Si-Si bonds break when increasing the thickness from 
30 to 100 /im. 

This phenomenon is even more visible in the inset of Figure [7] which shows the relationship 
{t^ / R) versus tp. The first two points which correspond to the 20 and 30 /im thick layers are on 
a straight line which passes by according to eq 14. Then, once the Si-Si bonds begin to break, 
the corresponding experimental points deviate strongly from this linearity. The results shown in 
Figures [H] and [7] are quite coherent. The shift of the adsorption isotherms (Figure O appears to 
have a bearing on the concomitant decreasing of the radius of curvature (Figure [7]). 

Whereas the 100 /im thick layer behaves like a membrane in adsorption-desorption process, 
it is not completely detached from the substrate, since the radius of curvature has a finite values 
(~ 1 m). To check whether these breakings are visible, we have performed scanning electron 
microscopy observations. The image of the porous layer-substrate interface shown in Figure [8] 
did not reveal any fracture at a scale of ~ 100 nm. 

When the constraint of the substrate is partially relaxed by the breaking of Si-Si bonds, the 
isotherms coresponding to the 100 /im thick supported layers are superimposed to that of the 
membranes. This indicates that the PSDs for both types of layers are similar, which justifies the 
assumption we made at the end of Section [3T] This suggests also that the difference between the 
two condensation branches observed for the thin layers (see Figure[T|) is not due to pore end effect 
but has the same cause as the difference between the adsorbed amounts before condensation. The 
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Figure 8: Scanning electron microscopy image in cross section of the porous silicon-silicon sub- 
strate interface for a porous layer 100 /im thick with a porosity of 50%. The amplitude of the 
undulations of the porous layer-substrate interface is about 100 nm, that is, 10^^ in relative 
terms. 



gap between the two condensation branches is entirely due to the different elastic deformations 
undergone by the thin supported layer and the membrane. 

Before discussing our results, let us comment recent experiments made on porous silicon layers 
similar to ours. 33J In this work, the adsorption isotherms for 30 fim thick supported layer and 
membrane have been found to be superimposed, which is in contradiction with our results. It is 
surprising that the authors who refer to our previous work published in 2002, [16j which shows, 
for the first time, that the isotherms for 20 fim thick supported layer and membrane are distinct, 
did not mention at all the disagreement between the two experimental results. 

How to explain such a disagreement? 

1) The process used to detach the porous layer from the substrate is not described in ref [33] 
but we suppose that it is the same as that currently used in the literature and by us, that is, an 
increase of the current intensity after the pores reach the desired length. We have seen above 
that this process lets the porous layer unchanged. Thus, in both experiments, the supported 
layer and membrane are identical. 

2) The (100) Si substrates used in both experiments are similar, 2 — 5 mil. cm in ref ^3] and 
2.7 — 3.3 mfi.cm in the present paper, like as the current density (20 mA/cm^). However, the 
etching solutions were different. In ref '33], the etching solution was HF(48%)/EtOH=l:l while 
in our case it was HF(40%)/EtOH=:3:l. The HF and water concentrations were hence different 
in the two experiments and consequently the porosity and etching rate which depend importantly 
on these parameters (see below). 

(i) Concerning the porosity, according to the calibrations we made, the conditions chosen 
in ref [33] should lead to porous layers with a porosity around 60%. Thus, the layers studied 
in ref 133] do not seem to be very different from ours which is coherent with the fact that the 
superimposed hysteresis loops for the membrane and supported layer are observed at pressures 
only slightly higher than represented in Figure 

(ii) The thickness of the porous layers was stated to be 30 fim but the authors did not tell 
us if they measured the thickness by weighing or if they chose the formation time from our 
calibrations. |15j The etching rate depends importantly on the water and HF concentrations. 
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For a given current density, for instance 50 mA/cm^, the etching rate for HF(40%)/EtOH=l:l 
is 1.7/im/min whereas for HF(48%)/EtOH=l:l it equals 2.5/zm/niin, that is 1.5 times higher. 
We have shown in this paper that the film thickness is a determinant parameter for breaking of 
Si-Si bonds like the stiffness and the effective lattice mismatch which depend on the porosity: 
for the porous layers we studied, the critical thickness from which the Si-Si bonds begin to break 
equals 30 fim. What is the critical thickness for the layers of ref [33]? In the absence of more 
informations we cannot say anything more than: if indeed the isotherms for supported porous 
layer and membrane in ref |33j are superimposed, the apparent disagreement between this result 
and ours represented in Figure [1] should have the same explanation as the apparent disagreement 
between our Figures [1] and [21 that is, the partial relaxation of the constraint of the substrate 
when making thicker layers. 

4 Discussion 

We have seen in this paper that the isotherms corresponding to porous layers supported by 
the substrate and the corresponding membranes are distinct whereas these layers are strictly 
similar. If we take into consideration the elastic deformation of the porous solid, these isotherms 
correspond in fact to different experiments. 

The porous silicon layer, when attached to the substrate, is submitted, during the adsorption 
process, to unilateral deformation perpendicular to the substrate, that is, along the pore axis 
(longitudinal deformations). The problem is equivalent to the classical case of the compression 
of a block held by lateral walls in such a way that its tranversal dimensions may be regarded 
as constant. When the porous layer is detached from the substrate to obtain a membrane, 
deformations perpendicular to the pore axis (transversal deformations) can also occur so that the 
extension along the pore axis is lower. The longitudinal deformations depend on the "effective" 
elastic constants, for the membrane and k^i for the supported layer, which connect the 
variation of the surface free energy with the deformations in the [100] direction. These elastic 
constants are related one to the other by the relation = Pkgi ~ 1.2ksi. The fact that /? ~ 1.2 
is close to 1, or, in an equivalent way, that Poisson's coefficient has a low value {vp ~ 0.1), 
reflects the high anisotropy of the porous system. Indeed, /3 = 1 and Vp = would correspond 
to a porous system whose transversal dimensions are unaffected by the deformation along the 
pore axis. On the other hand, for isotropic porous material such as porous glass, [28| taking 
j/p ~ 0.15 - 0.2 we obtain /3 > 2. 

The relaxation of the constraint of the substrate allows the porous system to jump from a 
state corresponding to onc-dimcnsional deformations to another state corresponding to tridimen- 
sional deformations. Although the transversal deformations are low compared to the longitudinal 
deformations since the two type of layers have very similar " effective" elastic constants along the 
[100] direction, the change of state corresponds to a significant decrease of the total free energy 

10%). This decrease results in a decrease of the adsorbed amount before condensation (10%) 
accompanied by a shift of the condensation pressure toward higher pressure, typically 0.7 — s- 0.73. 

Let us now consider a hypothetical porous system which could jump from an initial state 
corresponding to unilateral deformation as it is the case for the supported layer, to a final state 
corresponding to uniform tridimensional deformation as it is the case for porous glass or silica 
aerogel (/3 > 2). Our results suggest that, for the isotropic system, the decrease of the adsorbed 
amount before condensation should be much higher than 10% and the condensation pressure 
much higher than 0.73. 

In the light of this interpretation, let us analyze some results concerning other porous systems 
such as MCM-41 or SBA-15 silica, which are periodic arrangements of cylindrical pores. For 
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porous silicon and the SBA-15 material that we have studied, [13] the relative condensation 
pressures are similar, p ~ 0.7, while the mean pore size is ~ 13 nm for porous silicon and 
~ 8.5 nm for SBA-15 silica. In a previous paper, [TB] we had noted that the mean condensation 
pressure for porous silicon is lower than that expected from the classical Kelvin equation (inert 
adsorbent), roughly O.TPq instead of 0.9Po- Moreover, calculations based on density functional 
theory [331 (33 failed to reproduce such low condensation pressures. The condensation pressure 
for porous silicon seems hence "abnormally" low compared to that of SBA-15 silica having similar 
pore size and to theoretical values found under the assumption of an inert adsorbent. Note that 
to quantitatively reproduce the actual size for large pores in MCM-41 (6 — 7 nm) it is necessary 
to apply the Kelvin equation for a hemispherical meniscus to the capillary condensation branch. 
|36j This means that, as for porous silicon, the condensation pressure is lower than that expected 
from the classical Kelvin equation, i.e., inert adsorbent and cylindrical meniscus. 

What do we know about the elastic deformation of MCM-41 and SBA-15? For such porous 
materials, the only parameter which can be measured is the distance between the pore lattice 
planes, that is, the transversal deformations, so that we cannot compare the longitudinal and 
transversal deformations. Concerning the transversal deformations, there are still few results 
[3 [H] but they all suggest isotropic deformation. Moreover, Poisson's coefhcient values for these 
materials (TJ [37] are between 0.17 and 0.24, values close to that of homogeneous materials and 
importantly higher than that of porous silicon {vp ~ 0.1). We conclude that these materials 
react to the surface pressure exerted by the adsorbate in a similar way as isotropic materials. 
According to the above interpretation, this can partly explain the gap between the condensation 
pressures for porous silicon and SBA-15. 

Other parameters influence the condensation pressure such as the adsorbate-porous solid 
interactions or the pore morphology. We can note, for example, that the pores of porous silicon 
have a polygonal shape more or less regular whereas those of SBA-15 seem to be quite cylindrical. 
The PSD centered on ~ 13 nm in porous silicon layers under study corresponds to cylindrical 
pores having the same cross section surface area as the polygonal pores. If, as it is probably the 
case, the condensation in a given pore occurs when a liquid bridge is formed somewhere in the 
pore, the important parameter is the minimum distance between two walls. To estimate this 
distance, we have determined the PSD of the circles inscribed in the polygones: it is centered on 
11.5 nm. Thus, porous silicon is perhaps not equivalent to an assembly of cylindrical pores with 
a diameter mean value equal to 13 nm but rather to 11.5 nm. 

The theoritical calculations performed under the assumption of an inert adsorbent fail to 
reproduce the "abnormally" low condensation pressures for porous silicon, while the disagreement 
for materials submitted to isotropic strains (SBA-15) seems to be less important. This does 
not mean that the influence of elastic deformation is less important for material submitted 
to isotropic strains than for material submitted to anisotropic strains. Porous silicon (high 
anisotropic strains) and porous glass (isotropic strains) have very similar qualitative hysteretic 
behavior with similar boundary hysteresis loops and subloops inside the main loop. The difference 
between the total free energies for the condensation and evaporation branches for both systems is 
essentially due to the difference between the surface free energies, that is, the elastic deformations, 
which suggests that the height of the barrier to evaporation is essentially due to the elastic 
deformation. [3] For a given pore size, the position of the hysteresis loop on the pressure axis 
can depend on whether the porous system is submitted to anisotropic or isotropic strains, but 
the hysteretic behaviors are qualitatively similar. 

An alternative explanation based on the assumption of inert adsorbent has been proposed in 
the literature to account for the triangular shape of hysteresis loop in porous silicon, [ 33. .38, .39] 
while in SBA-15 or MCM-41 the shape is rather rectangular with steep and parallel condensation 
and evaporation branches. The hysteretic behavior of porous silicon would be due to the presence 
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of high disorder in the pores, fluctuation of pore diameter, [321 chemical heterogeneities, [35] and 
so on. One of the consequences of the hypothetical disorder introduced in each pore is that 
condensation is nucleated via the formation of liquid bridges so that pore ends do not play any 
role on the condensation process. Thus, the simulated isotherms for pores closed at one end (the 
supported layer) and for pores open at both ends (the membrane) are superimposed, [321 133] 
which is in contradiction with our experimental results. The disorder cannot explain why the 
experimental isotherms for the supported layer and membrane are distinct. 

On the other hand, we point out that our results suggest strongly that capillary condensation 
is nucleated via the formation of liquid bridges and not at the bottom of the pores. Indeed, 
we have seen that, in the region preceding capillary condensation, the normalized adsorbed 
amount at a given pressure decreases continuously with increasing layer thickness to reach the 
value corresponding to that of the membrane. We have explained this phenomenon by the 
partial relaxation of the constraint exerted by the substrate on the porous layer. Simultaneously, 
we observe that the condensation branch shifts toward higher pressures to finally superimpose 
itself on that of the membrane. The simultaneity of the two phenomena suggests that they are 
correlated. The question is whether the shift of the condensation branch observed when the 
thinnest supported layers are detached from the substrate (Figure could be partly due to pore 
end effect as modeled by Cohan. [T7] In the 100 /im thick supported layer, the major part of 
Si-Si bonds are not broken, since the curvature radius is small (1 m). Thus, the pores, or at 
least the major part of them, are closed at the porous layer-substrate interface. The pore size 
over pore length ratio is very small in the 20 /im thick layer, of the order of 5/10000, so that, 
if the formation of liquid bridges does not occur in the first 20 /im, there is no reason that it 
occurs at higher depth, since the average pore section does not vary with depth. Thus, the simple 
fact that the condensation branches for the 100 /im thick supported layer and membranes are 
superimposed suggests strongly that there is no pore end effect on the condensation process. 

The calculations of Naumov and co-workers, [33j reproduced a hysteresis loop with a triangu- 
lar shape similar to that experimentally observed in porous silicon. The geometry they considered 
consists of a slit pore divided in segments, with each segment having a diameter in the range 
4 — 8 lattice constant. For the porous silicon layers studied in the present paper, it would mean 
that each pore is composed of segments having diameter in the range 9—17 nm. In other words, 
the pore size distribution we determined from TEM images in plane view would represent the 
variation of the diameter of each pore. This should result in large undulations of the walls of 
each pores in opposite directions. These undulations should let the pore wall thickness remain 
unchanged, 5 — 6 nm, as well as the axes of the pores parallel one to the other. Indeed, we have 
shown by proton energy loss fluctuations measurements that the average direction of the pores 
is identical to within 0.1 of the [100] crystal axis. [40] We have never observed such undulations. 

The calculations presented in ref [331 are based on an unfounded hypothesis, and the fact 
that they reproduce the shape of the experimental sorption isotherms does not prove at all 
that the matter looks like the model, especially if these calculations are based on an unjustified 
assumption: an inert porous solid. 

To explain the presence of both a triangular hysteresis loop and an interaction mechanism 
in porous silicon where the pores are not interconnected, we follow a different way. We 
realized that the assumption of an inert adsorbent in adsorption process is unjustified. [^ Let 
us consider, for example, the variation of the total free energy at the top of the plateau between 
the condensation and desorption branches. We can see that the variation of the surface free 
energy, deduced from elastic deformation measurements, for porous glass and porous silicon, is 
of the same order of magnitude as the other terms of the total free energy and cannot thus be 
neglected. To neglect the elastic deformation, that is, the variation of the surface free energy, it 
would lead to a very small variation of the total free energy with regards to the chemical and 
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mechanical energies of the fluid. 



5 Conclusion 

The hnear deformation of a porous material during the adsorption-desorption process is pro- 
portional to the surface free energy, the proportionality factor depending on the elastic constants 
of the porous solid. [31 UHl E] This relation is the result of two equilibrium conditions, one 
thermodynamic and the other mechanic. [H] This means that a stress external to the porous 
layer can change the adsorbed amount and hence the condensation pressure. It is what we show 
in the present paper. The external stress is here exerted by the substrate on the porous layer. 
With porous silicon one can make identical porous layers either supported by or detached from 
the substrate. Whereas the two types of layers are identical, the isotherms are distinct. During 
gas adsorption in the supported porous layers, the planes perpendicular to the substrate are 
forced to have the same interatomic spacing as that of the substrate so that the porous layer 
is submitted to unilateral deformation, while in its relaxed state (the membrane) the porous 
system has three degrees of freedom and can find a new configuration for the solid atoms which 
significantly decreases the total free energy 10%) and the adsorbed amount (~ 10%) and 
increases the condensation pressure (0.7 0.73) compared to that of the supported layer. 

The stress exerted by the substrate can be partly relaxed by making thicker porous layers. 
The relaxed lattice parameter of a porous silicon membrane is greater than that of silicon, |20j 
so that the supported porous layer behaves like an epitaxial layer over silicon, which would have 
a positive lattice mismatch: the formation of a porous layer in a Si wafer produces an elastic 
bending of the wafer which then becomes convex. As the porous layer thickness is increased, 
the stress at the porous layer-substrate interface increases and the radius of curvature decreases 
although it remains higher than expected by the theory. This is due to the breaking of Si-Si 
bonds (dislocations) caused by the stress at the interface. The stress exerted by the substrate is 
thus partly relaxed by the formation of dislocations, which results in the shift of the isotherms 
toward that of the membrane when the porous layer thickness is increased from 30 to 100 fim. 

The longitudinal deformations (along the pore axis) depend on the "effective" elastic con- 
stants, km for the membrane and ksi for the supported layer, which connect the variation of 
the surface free energy with the deformations in the [100] direction. These elastic constants are 
related one to the other by the relation km = f3ksi ~ 1.2ksi. The large disagreement between 
this experimental value (/9 ~ 1.2) and that found in porous material submitted to isotropic de- 
formation {(3 > 2) indicates clearly that perpendicular to the pore axis the elastic deformation 
in the relaxed state of porous silicon (the membrane) is much lower than longitudinal deforma- 
tions. Although the transverval deformations are weak, the relaxation of the substrate results in 
a significant decrease of the total free energy and the adsorbed amount quoted above. 

This suggests that, for a porous system which would be subjected to an uniform tridimensional 
deformation {(3 > 2), the decrease of the adsorbed amount and the shift of the condensation 
pressure toward higher value should be much higher than those observed for the porous silicon 
membrane. This can partly explain why the condensation pressure for porous silicon layers is low 
compared to that of other ordered porous materials such as MCM-41 and SBA-15 silica having 
similar pore size, since these latter materials seem to react to the surface pressure exerted by the 
adsorbate in a similar way as isotropic materials. 

This does not mean that the influence of elastic deformation is less important for material 
submitted to isotropic strains than for material submitted to anisotropic strains. The hysteretic 
behaviors for porous silicon (anisotropic strains) and porous glass (isotropic strains) are very 
similar, and the height of the barrier to evaporation is essentially due to the elastic deformation 
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for both systems. The position of the hysteresis loop on the pressure axis depends on whether 
the porous system is submitted to anisotropic or isotropic strains, but the hysteretic behavior is 
quahtatively similar. 

The experimental results presented in this paper show the significant influence of stresses 
external to a porous system on the adsorption process. This supports recent Monte Carlo simu- 
lations [5] which show that the elastic strains affect the fluid's phase behavior. This supports also 
our thermodynamic approach and the interaction mechanism proposed, [5] in which the external 
mechanical stress is imposed on a given pore by the neighbors. 
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